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The transverse mass spectra of J /if) and ip' mesons and £1 hyperons pro- 
duced in central Au+Au collisions at RHIC energies are discussed within a 

. 

O . statistical model used successfully for the interpretation of the SPS results. 
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The comparison of the presented model with the future RHIC data should 
serve as a further crucial test of the hypothesis of statistical production 
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of charmonia at hadronization. Finally, in case of validity, the approach 



X 



should allow to estimate the mean transverse flow velocity at the quark 
gluon plasma hadronization. 

The concepts of chemical (hadron multiplicities) and kinetic (hadron momentum spec- 
tra) freeze-outs were introduced to interpret data on hadron production in relativistic 
nucleus-nucleus (A+A) collisions. The equilibrium hadron gas (HG) model describes re- 
markably well the light hadron multiplicities measured in A+A collisions at the SPS |1| 
and RHIC energies, where the creation of a quark-gluon plasma (QGP) is expected. 
Recently it was found that also charmonium yield systematics in nuclear collisions at 
the SPS @] follow the pattern predicted by the hadron gas model. The hadronization 
temperature parameter extracted from the fits to the hadron multiplicities is similar for 



both energies: T# = 170 ± 10 MeV. It is close to an estimate of the temperature Tc for 
the QGP HG transition at zero baryonic density. 

Experimental results on inclusive hadron spectra and correlations show evidence for 
a hydrodynamic expansion of the matter created in heavy ion collisions. Strong trans- 
verse flow effects in Pb+Pb collisions at the SPS (average collective transverse velocity 
is approximately 0.5 -j- 0.6) are firmly established from the combined analysis [|5]|| of the 
pion transverse mass spectra and correlations. The kinetic ('thermal') freeze-out of pions 
and nucleons seems to occur at a rather late stage of an A+A reaction. The thermal 
freeze-out temperature parameter of pions measured by the NA49 Collaboration for 
central Pb+Pb collisions at the SPS is T f = 120 MeV. 

Further exploration of idea of the statistical J / ip production || led to the formulation 
of the hypothesis that the kinetic freeze-out of J ftp and ip' mesons takes place directly at 
hadronization || . This means that chemical and thermal freeze-outs occur simultaneously 
for those mesons and they, therefore, carry information on the flow velocity of strongly 
interacting matter just after the transition to the HG phase. A possible influence of the 
effect of rescattering in the hadronic phase on the transverse momentum (p?) spectra 
was recently studied within a "hydro + cascade" approach fl,|10||. A+A collisions are 
considered there to split into three stages: hydrodynamic QGP expansion ("hydro"), 
transition from QGP to HG ("switching") and the stage of hadronic rescattering and 
resonance decays ( "cascade" ) . The switching from hydro to cascade takes place at T = T c , 
where the spectrum of hadrons leaving the surface of the QGP-HG transition is taken 
as an input for the subsequent cascade calculations. The results PJIH suggest that the 
Pt spectrum of Qs is only weakly affected during the cascade stage. The corresponding 
calculation for charmonia are not yet performed within this model, but a similar result 
may be expected due to their very high masses and low interaction cross sections. 



In previous work [11] devoted to the analysis of the SPS data it was demonstrated 
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that the measured transverse mass (mr = Jp\ + m 2 ) spectra of Or 11 hyperons |12| and 



J /ip and if)' mesons []13| produced in Pb+Pb at 158 A-GeV collisions can be reproduced 
within a hydrodynamical approach using the same freeze-out parameters: hadronization 
temperature T = 170 MeV and the mean transverse flow velocity Ut — 0.2. Within our 
approach the value of the vt parameter extracted in this way should be interpreted as 
the mean flow velocity of the hadronizing QGP. 

In the present letter we discuss the transverse mass spectra of J / ip and ip' mesons and 
Q hyperons in central Au+Au collisions at RHIC energies within the statistical approach 
to charmonium production successfully used for the interpretation of the SPS data |3|JTT|]. 

Assuming kinetic freeze-out of matter at constant temperature T, the transverse mass 
spectrum in cylindrically symmetric and longitudinally boost invariant fluid expansion can 
be approximated as [fLlp : 

dNi r R ( ttit cosh yj- \ /p^sinh|/r\ 
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where dt = tanh -1 vt is the transverse fluid rapidity, R is the transverse system size, 
Ki and Iq are the modified Bessel functions. The spectrum (0) is obtained under the 



assumption that the freeze-out occurs at constant longitudinal proper time r = \/t 2 — z 2 
(t is the time and z is the longitudinal coordinate), i.e. the freeze-out time t is independent 
of the transverse coordinate r. The analysis of the numerical calculations done within a 
hydrodynamical model |10] shows that the latter is approximately fulfilled. The quality of 
the approximation made gets better for heavy particles considered here because a possible 
deviation from Eq.(|l]) decreases with increasing particle mass at constant pr- In order to 
calculate ([!]) the function Dt{t) has to be given. A linear flow profile, yx{r) = y^ ax ■ r/R, 



is often assumed in phenomenological fits |14| and we used it as well. Recent numerical 



hydrodynamical calculations [IIJ show that deviations from this simplified profile are 



smaller than 10% for both SPS and RHIC energies. Their importance for our results is 



discussed later in the paper. 

In Figs. 1-3 the measured m T -spectra of Cl~ Hl2| , J/ip and ip' EH produced in 
Pb+Pb collisions at 158 A-GeV are shown together with the fit performed using Eq. ([!]). 



The model parameters are: T = 170 MeV, y™ ax = 0.28 [II]. The value of y™ ax found 
from fitting the SPS data appears to be close to the numerical estimate done within a 
hydrodynamical approach (yjp ax w 0.3) [jnj. This encourages us to use a corresponding 



estimate of y™ ax = 0.55 made for Au+Au collisions at ^ NN = 130 GeV [10] for the first 
guess of the mr-spectra of charmonia at RHIC energies within our approach. As follows 
from the analysis of hadron yields within the HG model, in our calculations we use 
a hadronization temperature T = 170 MeV, the same as in the case of the SPS data. 
The resulting m-p-spectra for Q, J/tp and tp' are shown in Figs. 1, 2 and 3, respectively. 
The normalization of the SPS and RHIC spectra is arbitrary. Due to large transverse 
flow velocity of hadronizing QGP at RHIC energies, the particle spectra are expected to 
significantly deviate from the simple exponential form of the m-p-distribution: 

dN 



C exp (- pi) , (2) 



rriTdmjn 

where C and T* are normalization and inverse slope parameters, respectively. 

In order to quantify the deviations of the RHIC spectra from Eq. (2) we plotted in 
Fig. 4 the local inverse slope defined as: 
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versus % — m for J/ip, ip' and Q. The T* parameter was calculated using Eq.(|l]) with 
the values of parameters used previously for the SPS and RHIC energies. No significant 
dependence of T* is observed in the SPS case y^ ax = 0.28, whereas a very strong decrease 
of T* with increasing in the region — m<l GeV is seen for J / ip and ip' at RHIC 
^yjnax _ q_55) anc } moderate one is expected for Q hyperons. Thus the shape of the 
m-p-spectrum is strongly dependent on the magnitude of the transverse flow and the mass 



of the particle. These characteristic "hydrodynamical" features of the spectrum should 
allow for a clear test of our approach in the near future. 



To be more specific, we follow the suggestion of Ref. |1(J and consider two domains 
< ttit — m < 0.6 GeV and 0.6 < my — m < 1.6 GeV of m^-spectra. Within these 
rriT domains the Q, J /if) and if)' spectra at RHIC calculated according to Eq. (1) were 
approximated by Eq. (2) and the values of C and T* were found by the maximum 
likelihood method. The inverse slopes in low ttit domain (long dashed lines in Figs. 1 - 
3) are larger than the high rriT ones (dashed-dotted lines in Figs. 1-3). However, the 
difference between low and high ttit inverse slopes strongly depends on the particle mass 
- for Q the inverse slopes differ by 13 per cent only, whereas for if)' by about 40 per cent. 

It is important to check the sensitivity of the above results to the assumed flow profiles 
(density and velocity) of hadronic matter at freeze-out. We first note that because of the 
basic assumption of our model (Q, J /if> and if>' freeze-out at QGP hadronization) the 
freeze-out density at the boundary between the mixed and hadronic phases is determined 
by the hadronization temperature T = 170 MeV and hence the density profiles of these 
particles are constant as the function of r. Thus the uncertainty is due to the velocity 
profile only. The numerical hydrodynamical calculations |H| indicate that the transverse 
rapidity profile may deviate from the linear dependence assumed above by less than 10%. 
In order to estimate a possible influence of these deviations we calculated the spectra of 
Q, J /if) and if)' using a profile of the form Ut{^) = y™ ax ■ (r/R) a for a = 0.9 and 1.1. In the 
calculations we keep the average transverse velocity as for y™ ax = 0.55 and a = 1 case. 
It turns out that the shape of the spectra is almost insensitive to the assumed value of 
a parameter in the whole range of transverse mass for Q and for ttit — m > 0.6 GeV for 
J /if) and ip' mesons. In the low my domain (my — m < 0.6 GeV) the fitted inverse slope 
parameters for J /if) and if)' mesons vary by 3 per cent and 5 per cent, respectively, under 
considered variations of the a parameter. Therefore we conclude - the principal results of 
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our analytical model are weakly affected by the approximations made. 

The transverse flow velocity of hadronizing matter depends on the energy and central- 
ity of the collision. Thus, in general, we consider the parameter y™ ax as a free parameter, 
who's value should be extracted from the fit of Eq.(|l]) to the measured spectrum. Provid- 
ing the fit is successful, the fitted yjp ax may be treated as an estimate of the transverse 
flow velocity of hadronizing matter. 

The rriT spectra in Au+Au collisions at RHIC were also calculated within a hydrody- 



namical ansatz in Ref. [15|. There is however a principal difference between our approach 



1 1] and the one used in Ref. |I5|| . We consider the early kinetic freeze-out for f2, J/ip 
and ip' only, whereas in |TS| the simultaneous (chemical and thermal) freeze-out was pos- 
tulated for all hadrons. This generates a substantial difference in the model predictions 
for the RHIC freeze-out parameters. 

The results on charmonium production in central Au+Au collisions at RHIC energies 
are expected to be available very soon. They will, obviously, allow to test various ap- 
proaches |16H20| developed in the course of the analysis of the corresponding SPS data. 
In particular, a statistical approach to charmonium production in nuclear collisions yields 
well defined and characteristic predictions for charmonia multiplicities || and the trans- 
verse mass spectra discussed in this letter. The analysis of the m^-spectra of J/ip and 
ip' mesons and Q hyperons within our approach should yield a unique information on the 
transverse flow velocity of hadronizing QGP in A+A collisions at RHIC energies. 
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FIG. 1. The transverse mass spectra (in arbitrary units) of 0~ hyperons produced in 
Pb+Pb (Au+Au) collisions. The points indicate experimental results measured for Pb+Pb 
collisions at 158 A-GeV. The dashed line shows the fit of the hydro dynamical model to these 
data (T = 170 MeV and y™ ax = 0.28). The prediction of our approach for central Au+Au 
collisions at *fs NN = 130 GeV is indicated by the solid line (T = 170 MeV and yf ax = 0.55). 
The long dashed and dashed-dotted lines correspond to the exponential spectrum fitted in the 
intervals (mr — m) £ [0, 0.6] GeV and (my — m) £ [0.6, 1.6] GeV), respectively. 



9 



10 



10 



10 



10 



10 



-10 



10 



-12 







T* = 589 MeV 
T* = 399 MeV 




RHIC 



m T -m (GeV/c ) 



FIG. 2. The same as in Fig. |T| but for J/tp mesons. 
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FIG. 3. The same as in Fig. [l| but for ip' mesons. 
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FIG. 4. The dependence of the inverse slope parameter T* of J /ip (solid line) and ^ (long 
dashed line) mesons and hyperons (short dashed line) on tut — m expected within the model 
of statistical production at QGP hadronization. The left panel indicates the behaviour given by 
the fit of the hydrodynamical ansatz to the data on Pb+Pb collisions at 158 A-GeV, whereas 
the right panel shows the predictions for central Au+Au collisions at v^TViv = 

130 GeV. 
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